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STRESS -INTENSITY FACTORS BY BOUNDARY COLLOCATION FOR SINGLE- 


EDGE-NOTCH SPECIMENS SUBJECT TO SPLITTING FORCES 
by Bernard Gross and John E. Srawley 
Lewis Research Center 


SUMMARY 

A boundary -value-collocation procedure was applied in conjunction with the 
Williams stress function to determine values of the stress -intensity factor Kj 
for single-edge cracks in plate specimens subject to splitting forces applied 
close to the crack and acting transversely to it. The results are presented in 
terms of the dimensionless quantity Y = KjBH^/ 2 /pa^ where B and H are the 
specimen thickness and half -depth, a is the effective crack length, and P is 
the applied load. The results are practically independent of the ratio of ef- 
fective specimen length to specimen half -depth q/H, when this ratio is not less 
than a/H + 2, and are then In excellent agreement with those derived by other 
investigators from compliance measurements. In the limit, as H/a approaches 
zero, the value of Y approaches that obtained in an elementary analysis which 
treats the specimen as a pair of built-in cantilever beams. 


INTRODUCTION 

Plate specimens of the general configuration shown in figure 1 are used for 
plane strain fracture toughness measurements (ref. l) and for fatigue crack 
propagation rate measurements . This type of specimen is essentially a narrow 
plate with a crack extending from one of the short edges along the longitudinal 
centerline. The forces P are applied close to the crack and act transversely 
to it, thus tending to split the specimen by extending the crack. Such speci- 
mens are referred to in this report as single -edge -notch specimens subject to 
splitting forces. They are sometimes referred to elsewhere as crackline loaded 
specimens, or as double cantilever specimens. The length of the crack a is 
taken as the distance from the end of the crack to the line of action of the 
load P. (All symbols used in this report are defined in the appendix.) 

A first approximation to the strain energy release rate with crack exten- 
sion for such a specimen is readily calculated from simple beam theory by re- 
garding the specimen as a pair of built-in cantilever beams. This approach was 
used by Gilman in a study of the mechanics of cleavage (ref. 2). In terms of 
the stress -intensity factor Kj, the result obtained is that Kj is equal to 



2 V3 Pa /h 3 /2 for small values of H/a, where H is the depth of each canti- 
lever arm (ref. 3) . A better approximation was obtained by Ripling, Mostovoy, 
and Patrick (ref. 4), who modified the cantilever beam treatment and used ex- 
perimental compliance measurements to establish the best value for a correction 
term. 


Because of the potential usefulness of this type of specimen, it was con- 
sidered desirable to conduct an independent study of the configuration by the 
boundary-collocation procedure used in earlier studies of other types of single 
edge-notch specimens by Gross, Srawley, and Brown (refs. 5, 6, and 7). It was 
decided that this study should also include an evaluation of the effect of a 
variable not considered in reference 4, namely, the ratio of the effective spec 
imen length to the effective crack length q/a. It is apparent that the stress 
intensity factor will be independent of q/a for sufficiently large values of 
that ratio, but it is of practical interest to know how large a ratio is suffi- 
cient. 

In general, the stress -intensity factor Kj will depend on the five vari- 
ables P, a, q, H, and B (fig. l) . It is convenient to express the relation 
between Kj and four of these variables in dimensionless form in terms of the 

stress -intensity coefficient Y = KjBH 3 /^/pa, which is a function of the ratios 
a/H and q/H. For practical purposes, the dependence of Y on q/H is neg- 
ligible when q/H is sufficiently larger than a/H. The relation between Y 
and a/H for sufficiently large q/H will be referred to as the K-calibration 
for the single-edge-notch specimen subject to splitting forces. The K- 
calibration results discussed in this report cover a range of values of a/H 
from 2 to 10, which is sufficient for practical purposes. For each value of 
a/H, results were computed for values of the parameter q/H of a/H + 1, 
a/H + 2, and a/H + 3 sufficient to define the practical limit of influence of 
this parameter. For a/H equal to 4, results were also computed for q/H 
equal to a/H + 6. 


ANALYTICAL AND COMPUTATIONAL PROCEDURE 


The method of analysis consists in finding a stress function X that sat- 
isfies the biharmonic equation V^X =* 0 and also the boundary conditions at a 
finite number of stations along the boundary of a single-edge-notched specimen, 
such as shown in figure 2. The biharmonic equation and the boundary conditions 
along the crack are satisfied by the Williams stress function (refs. 8 and 9). 
Because of symmetry (fig. l), the coefficient of the sine terms in the general 
stress function must be zero, hence 
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The stresses in terms of X obtained by partial differentiation are as follows: 
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The boundary -collocation procedure (ref. 5) consists in solving 2m si- 
multaneous algebraic equations corresponding to the known values of X and ei- 
ther 5x/chc or 5x/By at m selected stations along the boundary ABCD of 
figure 2; thus, values for the first 2m coefficients of the Williams stress 
function are obtained when the remaining terms are neglected. Only the value 
of the first coefficient d^ is used for the present purpose since the stress - 

intensity factor is equal to -•^2-rt dj, as shown in reference 5. 

The required values of X and its first derivatives at the m selected 
boundary stations were obtained from distributions of bending moment and para- 
bolic shear stress distribution (fig. 2) equivalent to the actual loads shown 
in figure 1. The equations for these boundary values in dimensionless form 
are as follows : 
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The distance V of the end boundary from the crack tip (fig. 2) was chosen 
to be 0.8 a. From physical considerations, it is clear that the boundary should 
be chosen neither close to the crack nor close to the load point because of the 
stress-field disturbances near these positions. Preliminary studies similar to 
those described in reference 5 established that V/a equal to 0.8 was about 
optimum, and that minor variations from 0.8 had a negligible effect on the K- 
calibration. 

For each set of selected values of the primary variable a/H and the 
parameter q/H, the collocation computation was carried successively over the 
range of m equal to 21 to 48 stations in increasing increments of three. Fig- 
ure 3 is fairly representative of the resultant variation of the K-calibration 
against the number of boundary stations. Except for one case, the variation of 
the dimensionless stress -intensity factor Y obtained for m equal to 42, 45, 
and 48 was less than 1 percent and in most cases it was a small fraction of a 
percent. For a/H equal to 2 and q/a equal to 5/2, the variation was approx- 
imately 1^- percent. From this observation taken together with the nature of the 

trends of the Y values, the K-calibration results shown in table I were based 
on the average value of Y for m equal to 42, 45, and 48. 

RESULTS AID DISCUSSION 

Table I shows a comparison of results of the present work with correspond- 
ing results from reference 4 in terms of the dimensionless quantity KjBH 3 / 2 /Pa 

as a function of relative crack length a/H. Over the range of relative crack 
lengths a/H varying from 2 to 10, the effect of relative specimen length q/H 
becomes negligible when not less than a/H + 2. Figure 3 shows a typical result 
in evaluating the effects of a large q/H ratio, a/H + 6, on the dimensionless 
quantity Y. The agreement between the analytical results for q/H greater 
than or equal to a/H + 2 and the experimental results of reference 4 is within 
1.3 percent for all cases. 

For a/H less than 2, from physical considerations one can conclude that 
an interaction between the stress field of the applied load and the crack tip is 
likely to occur. Thus, use of simple beam theory to obtain the stress function 
and its associated derivative at location V, would be an oversimplification of 
a complex condition. 

The plot of Y against the reciprocal of relative crack length (fig. 4) is 
well fitted by a straight line passing through Y = 3.46 at H/a = 0. The 

value 3.46 is equal to the numerical coefficient 2-/3 obtained as a first 
approximation from simple beam theory, as discussed in the introduction. Thus, 
extrapolation of the collocation results to the limiting case of a very long, 
slender specimen leads to a value which is in good agreement with the elementary 
treatment. 
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MEASUREMENT CAPACITY IN RELATION TO SPECIMEN DEPTH 


The currently accepted criterion for evaluating Kj c , the critical value 

of K at point of instability of crack extension in first or open mode, is 
that the nominal stress at the crack tip should not exceed the yield strength of 
the material Oj S in a valid Kj c test. Thus, the maximum nominal stress 

6Pa/BH^ should not exceed cr ys , which, on substituting in expression 

Y 2 = K 2 B 2 H 3 /P 2 a 2 gives = Y 2 /72, where Cjk is used to denote the 

maximum acceptable value Kj c . Figure 5 shows the dependence of dimensionless 
Kj c measurement capacity on the relative crack length a/H for a sufficiently 
large effective specimen length. The specimen measurement capacity decreases 
with increasing relative crack length. For relative crack lengths greater than 
6, the drop in capacity becomes exceedingly small. 

The dependency of Kj c measurement capacity on the specimen thickness is 
discussed at some length in reference 1 and is unrelated to the present work. 


Lewis Research Center, 

National Aeronautics and Space 
Cleveland, Ohio, December 


Administration, 
9, 1965. 
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APPEKDIX - SYMBOLS 


a 


B 


C 


IK 


d 2n' d 2n-l 

H 

K I 


m 

P 

q. 

r 

V 


Y 

G 


'x 


ys 


xy 


x 


effective crack length 
specimen thickness 

estimate of maximum value of K^ c that can he measured with speci- 
men of given dimensions and yield strength 

coefficients of Wi ll iams stress function 
depth of arm of split beam 

stress -intensity factor of elastic stress field in vicinity of crack 
tip 

critical value of K at point of instability of crack extension in 
first or open mode, a measure of plane strain crack toughness of 
material 

number of selected boundary stations for collocation solution 

load applied to both ends of split beam 

effective specimen length 

polar coordinate referred to crack tip 

distance from crack tip to boundary for collocation analysis 
Cartesian coordinates referred to crack tip 
dimensionless stress -intensity coefficient, KBH 3 / 2 /pa 
polar coordinate referred to crack tip 
stress component in x-direction 

stress component in y-direction 
0.2 percent offset tensile yield strength 
shearing stress component 
stress function 
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TABLE I. - COMPARISON OF RESULTS OF PRESENT WORK 
WITH CORRESPONDING RESULTS FROM REFERENCE 4 IN 
TERMS OF THE DIMENSIONLESS QUANTITY KjBH 3 / 2 /^ 


AS FUNCTION OF RELATIVE CRACK LENGTH 



P 


Figure 1. - Edge-cracked plate specimen subject to splitting forces. 
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Number of boundary-collocation stations, m 


Figure 3. - Value of dimensionless stress-intensity coeffi- 
cient Y against number of boundary-collocation stations 
for relative crack length a/H - 4. 
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Reciprocal of relative crack length, H/a 

Figure 4. - Dimensionless stress-intensity coefficient as 


function of reciprocal of relative crack length. 



Relative crack length, a/H 

Figure 5. - Dependence of K Ic measurement 
capacity C IK of single-edge-notched speci- 
mens subjected to splitting forces on relative 
crack length. 
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